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1. Introduction 



Thermodynamics, or thermodynamic method, has been invented as a way to describe the 
exchange processes of energy and matter occurring at the molecular and atomic levels without 
considering the details of molecular motion. The very beginning of it has been inspired by 
steam engines efficiency and the initial context of thermodynamics concerned the macroscopic 
systems in equilibrium and/or quasi-equilibrium states. Statistical physics has been appeared 
as an attempt to describes the same processes with taking into account atomic hypothesis and 
atomic dynamic parameters, i.e. the velocities and the positions of atoms. The probabilistic 
language appeared to be very suitable for these purposes. But again the studying has been 
restricted by the macroscopic systems in quasi-equilibrium states mostly because of using 
variational principles like the maximum entropy principle. 

Today the relevance of thermodynamic formalism and applicability of statistical physics are 
questionable when the nanosystems, the systems far from equilibrium and the systems with 
strong interactions begin to be studied. First of all the problems appear for the definition of 
the temperature that is the key concept in the formalism of thermodynamics and statistical 
physics. The assumptions used to define the temperature have to be treated very carefully in 
the cases of nanosystems, systems with strong interactions and other complex systems 1 . 
The goal of the present paper is to remind the main conditions which has to be satisfied to 
introduce the physical quantity "temperature" and to discuss the possibilities for introducing 
the temperature in complex systems. As an example the model system - the system of 
interacting spins at external magnetic field - will be used to demonstrate the advantages and 
restrictions of using "spin temperature" concept. 



2. What is the temperature? 

The thermodynamic or statistical definition of the temperature can be found in any standard 
physics textbook. Here we reproduce briefly these procedures with emphasizing some points 
which are usually considered as the given. As far as the exchange processes of energy and 

1 During the preparation of the present paper the very detailed book has been published iBiro, 2011) 
where the basic concepts at the very heart of statistical physics are presented and their challenges in 
high energy physics are discussed 
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matter are the subject of investigations two macroscopic systems in equilibrium are usually 
considered. If one allows only energy exchange between these two systems the equilibrium 
means the equality of some physical parameter in this case. We can call this parameter as a 
temperature. But how we can measure and/or calculate this parameter? It is necessary to 
mention here that the key-point in the definition of the temperature is the existence of the 
so-called thermal equilibrium between two systems. 

In the framework the phenomenological approach - thermodynamics - the temperature is 
measured by the monitoring other physical parameters (expansion coefficient, resistivity 
voltage, capacity etc.). That is why there are so many different kinds of thermometers. The 
procedure of temperature measurements consists of the thermal contact (energy exchange) 
between the system under consideration and a thermometric body the physical state of which 
is monitored. This thermometric body should be as small as possible in order to do not 
disturb the state of the system during measurement. In the realm of very small systems such 
a procedure is rather questionable. What the size should be for the thermometer to measure, 
for example, the temperature of a nanosystem? Should the thermometric body be an atom or 
elementary particle in this case? But the states of atoms and elementary particles are essential 
quantum ones and can not be changed continuously. The excellent treatment of the more 
sophisticated measurements of temperatur e (spectral t emperature and radiation temperature) 
the reader can find in the very recent book ^Birol l20lTb . 

The simplest way to determine the temperature in statistical physics is based on the 
consideration of possible microscopic states for the given macrostate which is determined 
by the energy of system or other constraints. The number of possible microscopic states 
O(E) is extremely fast rising function of the energy of system E. For two systems being in 
the thermal contact the total number of possible microscopic states is given usually by the 
following product: 

n 12 = n(Ej) • n(E 2 ). (i) 

An implicit assumption in the notion of thermal contact is that the system-system interaction 
is vanishingly small, so that the total energy E is simply given by 

E = Ei + E 2 = const. (2) 

The product n(Ej) • Cl(E — Ei) shows the very sharp maximum (see Figured} and it is more 
convenient to study the extremal conditions for the logarithm ln{Cl{E\) ■ 0(E — Ej)) from 
which on immediately get the definition of the statistical temperature: 

M"("(Ei)) = dln(Cl(E 2 )) =R = ±_ m 
dEi dE 2 P k B T' { ' 

where kg is the Boltzmann's constant and T denotes the absolute thermodynamics 
temperature. In fact the taking of logarithm leads us to some additive quantity, and it is 
the property which is carried by the Boltzmann definition of entropy: 

S = k B lnCl{E). (4) 

Considering the system in the contact with the thermal bath (thermal reservoir) the same 
assumption about neglecting system-bath interaction leads to the existence of canonical 
(Gibbs) distribution for the probabilities to find the system in the state with energy E K : 



P(E ci ) oc exp(-pE a ). 



(5) 
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Fig. 1. Two systems in the thermal contact and the number of possible microscopic states for 
them and for the joint system in dependence on energy. Note that in the definition of thermal 
contact the energy of interaction between systems is vanishingly small, so the total energy is 
just E = Ei + E 2 - 

We have seen that all obtained results are valid only if one neglects by the system-system 
or system-bath interactions. In fact such neglecting proved to be the consequence of 
the so-called thermodynamic limit which is reached as the number of particles (atoms or 
molecules) in a system, N, approaches infinity. For the most systems in the thermodynamic 
limit the macroscopic extensive variables (energy, entropy, volume) possess the property of 
additivity like Equation[2] It is necessary to point out here that the possibility to describe the 
thermodynamic behavior of the system under consideration by the statistical physics methods 
as N tends to infinity is not automatically granted but depends crucially of the nature of 
the system. It has been shown many times that the Gibbs canonical ensemble is valid only 
for sufficiently short range interactions and there are examples - self-gravitating systems, 
unscreened Coulomb systems - for which the assumed additivity postulate is violated. 
Also it should be m entioned th at it is possible to give an alternative definition of the 
thermodynamic limit jBirol [201lh . If a system is so large that it itself can serve as a perfect 
thermal reservoir (bath) for its smaller parts, then one can consider this system as being in 
the thermodynamic limit. This definition is not restricted to large volume and large particle 
number. But again here the energy of spin-bath interactions is neglected and, in practice, such 
definition also implies the existence of additive quantities. 

The more rigorous definition of temperature in statistical physics is based on the maximization 
of entropy for the system composed from two subsystems (like we consider above) while 
the energy, volume, par ticle numbe r etc. are composed from the corresponding subsystems 
value (see, for example, ( Bir ot |201lJ) ) . Mathematically in denotes that one should look for the 
maximum of the entropy of the system 

S(E,V,N,...) =max (6) 
when the following equations are satisfied 



E = E-i @ E 2 
V = V 1 ffi v 2 
N = Ni © N 2 . 



(7) 
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Strictly speaking the composition law © in Equations [7] is not restricted only by addition. But 
finding the maximum of entropy[8] 

*S(E, V,N,. ..) = (g.^ + §- d E 2 ) + ( * + + (*Lm + §- 2 iN 2 ) + . . . = 

(8) 

when the total energy, volume and particle number are fixed, i.e. 

dE = (™ dEl + ll dE2 ) = o (9) 



y dEi dE 2 
f dV dV \ 

d y=(- d y 1 + - d y 2 )=o (io) 

(please note, that here energy, volume and particle number are considered as independent 
variables!!!) one gets the following expressions: 



as 1 \ /as i 



as i \ f ds i 



as i \ / as i 



(12) 
(13) 
(14) 



The left (right) part of Equation[T2lwill depends on the quantities describing subsystem 1 (2) 
only in the case, when the simplest addition is taken in[7]as the composition law. So again we 
conclude that the rigorous definition of the temperature can be done only for the extensive 
systems in thermodynamic limit. 

Now we will demonstrate that the thermodynamic limit does not exist for the systems with 
strong enough interactions and how this limit can be restored by some averaging procedure. 
We consider a N-particle system describing by the following Hamiltonian: 

N N 

u = Y J u m + Y J u lj , (is) 

;=i i<j 

where T-i^j is the Hamiltonian of free particle and Ujj describes the interaction between ;-th and 
/-th particles. It is obviously that because of second term in Equation[l5]the energy of system 
is not additive and the system can not be trivially divided into two- or more independent 
subsystems. By increasing the number of particles in the system in n times the energy of 
systems is not increased in n times too. The conditions of thermodynamic limit are violated 
here. But if the interaction between particles W,y is a short-range one and a rather small the 
concept of the mean field can be introduced, so the Hamiltonian|15lcan be re-written as: 
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N N 
i=l i=l 

So the energy becomes additive and the conditions of thermodynamic limit are satisfied. 
We note here that for nanoscale systems in which the contribution to their energy (or other 
quantity) from the surface atoms is comparable with that from the bulk volume atoms 
the non-additivity exists forever. The surface energy is increased in n 2 ^ 3 times when one 
increases the number of particles in n times. That is why the question about the definition 
of temperature for nanoscale system is very intriguing not only from the point of view of its 
physical measurements but a lso from the estimatio ns for the minimal length scale on which 
this intensive quantity exists iHartmann et al.l,l2004h . 

Last decades the non-extensive thermodynamics has being developed to describe the 
properties of systems where the thermodynamic limit conditions are violated. It is not 
a purpose of this paper to give one more review of non-extensive thermodynamics, its 
methods and for malism. The reader can found it in numerous papers, reviews and books 
(see, for example, lAbe & Okamotol l200ll; lAbe et all |2007|; iGell-Mann & Tsallisl |2004|; iTsallisl 
l2009f) ). Here we would like to underline only that the definition of temperature is very close 
related to the existence of thermal equilibrium and is very sensitive to the thermodynamic 
limit conditions, so a researcher should be very careful in the prescribing the meaning of 
temperature to a Lagrange multiplier when entropy maximum is looked for. 
In the next part of this work we will show how the thermodynamic formalism can be applied 
for the system of spins in an external magnetic field and will discuss the existence of two spin 
temperatures for one spin system. It is a real good example when one has to remember all 
conditions being used to 

3. Spin temperature 

An important progress in the description of the behavior of spin degrees of freedom in 
solids was reached with the use of thermodynamic approach. The relation between electron 
and nuclear magnetism on the one hand and thermodynam ics on the other hand was first 
established by Casimir and Du Pre jCasimir & Du Pri. Il938h . They introduced the concept 
of spin temperature for a system of non-interacting spins in an external magnetic field. But 
for a long time spin temperature was being considered an elegant theoretical representation 
only. Further investigations of spin thermodynamics showed that spin temperature can be 
experimentally measured and its change is connected with the transfer of heat and change 
of entropy. For example in the framework of spin temperature concept Bloembergen et al. 
iBloembergen et all Il948l; Il959t) built the classical theory of saturation and cross-relaxation 
in spin systems . The excellent account of spin tempera ture concept the reader can find in 

jAbragaml,l2004l:lAtsarkin & RodakLll972l;lGoldmanl,ll97(]h. 

In the classical theory of Bloembergen et al. iBloembergen et all Il948|) the Zeeman levels 
are considered as being infinitely sharp thus neglecting the broadening due to spin-spin 
interactions. This is only justified in cases of liquids and gases, where the rapid motion of the 
atoms or molecules averages the spin-spin interactions to zero and the spins can be considered 
as being independent of each other. However in solids the spin-spin interactions are normally 
so strong that the whole ensemble of spins acts as a collective system with many degrees of 
freedom. The next mai n step in the understanding of s pin thermodynamics in solids was 
made by Shaposhnikov jShaposlmikovl.ll947l;ll948l;ll949h . whose works were far in advance 
of the experimental possibilities of their verifying. In these little known works Shaposhnikov 
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pointed out the significance of taking into account the interactions inside spin systems. As 
long as the time taken for internal thermal equilibrium to be established is finite then to 
describe any state of spin system one needs to determine two thermodynamic coordinates, 
which are magnetization and spin temperature. As distinct from the theory of Casimir and 
Du Pre these coordinates are not connected between each other by Curie law. Moreover, the 
spin temperature and total energy of spin system were proven to be just thermodynamically 
conjugated variables, and the magnetization characterizes a state of spin system in an external 
magnetic field. 

Independently in 1955 to investigate the saturation in system of interacting spins Redfield 
(Redfield, 1955) introduced the hypothesis that under strong saturation the whole spin system 
stayed in internal equilibrium, thus permitting its description by one single temperature. 
The problem for an arbitrary degree of saturation was solved in 1961 by Provotorov 
who developed Shapos hnikov's id e as abo ut takin g into account spin- spin interactions. 
In well-known works dProvotorovl Il962alfbl Il963ft Provotorov showed that under some 
conditions an energy of spin-spin interactions that are small compared with the interaction 
with strong external magnetic field can be extracted into a separate thermodynamic subsystem 
called by the reservoir of spin-spin interactions. This thermodynamic subsystem has its 
own temperature different from the temperature of spin system in an external magnetic field 
determined in Gasimir and Du Pre theory. Therefore according to the Provotorov's theory any 
state of spin system can be described by two temperatures. The concept of two temperatures 
turned out to be fruitful and it was confirmed experimentally in electron as well as in nuclear 
magnetism. This concept led to the revision of some representations in the theory of magnetic 
resonance and relaxation and to the prediction of a number of unexpected physical effects. 
Today the two-temperature formalism being unusual from the point of view of statistical 
thermodynamics forms the basic framework for the theory of magnetic resonance in solids. 
But because of conceptual and mathematical difficulties the Provotorov's theory was well 
developed only in the so-called high-temperature approximation when the heat energy of 
spin greatly exceeds the energy of spin in external magnetic field. 

The attempts to extend the theory towards low temperatures (the energy of spin in external 
field is more than heat energy) demonstrated the principal difficulties in the choosing 
of thermodynamic variables and in the understa nding of energy redistribution inside 
spin-system (see, for review jTayurskiil Il989l: Il990h and references therein). But there are 
many experiments at low temperatures to interpret of which it is necessary to have a theory 
describing the spin thermodynami cs and kinetics in this case. Among them we point out 
dynamic polarization e xperiments jAbragam & Goldman!. 119821) . nuclear ordering in solids 
iBonamour et alj Jl99ll ): lLounasmaal jl989h . The object of this paper is to give a description 
of some theoretical approaches to the studying of spin-system in solids at low temperatures. 
In the next section we will summarize the concept of spin temperature, as far as it is very 
important for the understanding of low temperature thermodynamics of spin-system in 
solids. 

3.1 High and low temperatures 

If an electron (or nuclear) spin S is subjected to a magnetic field Hg in the direction z then the 
Zeeman interaction 

H z = to S z (17) 

establishes a set of 2S + 1 sublevels with energy E m = mcog (here and further the units fi = 1 
are used ); cvq = jHq is the?Larmor frequency, j is the gyromagnetic ratio, m is the magnetic 
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quantum number ( m = —S, —S + 1, . . . S — 1, S). In Equation [17] one must take the negative 
sign if S is a nuclear spin. For an ensemble of N identical spins S = 1/2 we can introduce the 
Zeeman level populations «+ and «_ ( n+ + n~ = N ), where w+ is the number of spins in the 
state m = +1/2 and n_ is the number of spins with m = —1/2. If the spins are in equilibrium 
with lattice the distribution of the different spins over the magnetic levels is determined by 
the Boltzmann law 



(18) 



where Tj_ is a lattice temperature. Such an equilibrium can be violated, for example, after 
saturation by radio-frequency field. In this case one can describe spin system introducing 
spin temperature T s that is distinguished from lattice temperature and is determined by the 
expression 



— = exp{— -} 



(19) 



It should be noted that in the case of thermal equilibrium with the lattice the levels of highest 
energy are always less populated than the lower ones. At the saturation it is possible to create 
situations in which the highest levels are more populated than the lower ones. If one assumes 
Equation[l9]to be still valid this situation corresponds to negative spin temperatures. The spin 
system is not then of course in equilibrium with lattice because negative temperature is only 
defined for systems with an upper bound in the energy spectrum. The Figure|2]illustrates the 
concept of spin temperature for the case S = 1/2 and positive, negative and infinite values of 
temperature. 
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Fig. 2. Illustration of spin temperature concept for the case S = 1/2. In the lower row the 
populations of level are shown for the case of positive, negative and infinite spin 
temperatures. 

In this part we shall discuss the thermodynamics of spin systems at low temperatures. We 
shall define "low temperature" as a temperature at which the energy of spin in strong external 
magnetic field is of the same order as the average heat energy or exceeds it, i.e. the following 
condition is true: 



<^0 > k B T L 



(20) 
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In the case of localized electron spins in insulators this condition denotes a magnetic field of 
about 50 kG and a temperature about 1 K. Concerning the energy of spin-spin interactions the 
assumption is made that it is small compared to the energy of spins in the external magnetic 
field. As we are not interested in the effects connected with phase transition into magnetically 
ordered state then we can assume that 

E ss « k B T L (21) 

where E ss is an average energy of spin-spin interactions. The temperatures at which an 
average heat energy of spins becomes more than the energy of a spin in an external field are 
called "high temperatures". The transition from the range of high temperatures into the range 
of low temperatures is accompanied first of all by the essential changes of the thermodynamic 
properties of spin system. 

3.2 High-temperature thermodynamics of spin system 

Let a consider the regular lattice of spins in an external magnetic field. The corresponding 
Hamiltonian of spin system has the following form 

Us = H Z + H ss (22) 

where 

«z = Wo£S? (23) 

i 

is the Hamiltonian describing the Zeeman interaction of spins with the constant magnetic 
field directed along z-axis, SJ is a longitudinal component of ;-th spin. In Equation [22] Hss 
describes the interactions of spins between themselves. For the sake of simplicity we will not 
take into account further interactions such as hyperfine interactions or interaction with crystal 
field. Temporarily we omit the interactions such as spin-phonon interactions and interactions 
with other external fields. The energy levels of the Hamiltonian H z are strongly degenerated 
because for any eigenvalue of the Hamiltonian one can find many combinations in which 
the eigenvalues of operators SJ can be taken. This degeneracy is removed by the spin-spin 
interactions. Further we will suppose the external field to be large compared to the internal 
fields induced by spin-spin interactions. Therefore the Hamiltonian of spin-spin interactions 
Hss can be considered as a perturbation. In the first order of the perturbation theory only 
those terms of H ss which don't cause the change of any eigenvalue of the Hamiltonian H z 
will give a contribution into the splitting of the energy level corresponding to this eigenvalue. 
So in the first order in the perturbation only part of Hss that commutates with H z will give 
a contribution in the broadening of the Zeeman levels. Usually this part is called the secular 
part and is presented as (4) 

Hss = \ E( A i) S ,^ + B ij S t S j ) ( 24 ) 
V 

where S ( + and are the transverse components of /-spin. Further we shall use the notation 
Hss exactly for the secular part of the Hamiltonian of spin-spin interactions. In this section we 
don't need the explicit form of spin-spin interaction constants A« and B,-y . Thus the energy 
spectrum of the Hamiltonian H s represents quasi-continued equidistant bands of energy 
levels and these bands are separated in the energy by equal intervals coq . 
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The distribution of spins over these levels is the result of spin-spin interactions during some 
characteristic time T2 that is of order of inverse width of the magnetic resonance line shape. 
As it follows from the experience usually in solids this time T2 is significantly shorter than 
spin-lattice relaxation time T-y. The distribution of spins over the energy levels entirely 
determines the thermodynamic state of a spin system. It was appeared to be very convenient 
to introduce two spin temperatures for describing this distribution. 

Formally the idea of two spin temperatures can appear from the following considerations. As 
long as spin system placed into strong external magnetic field has two first integrals - they 
are the Zeeman energy and the energy of the secular part of spin-spin interactions - then in 
order to describe a sufficiently slow evolution of any thermodynamic state it is natural to 
choose these first integrals as the thermodynamic coordinates. Such a choice becomes very 
attractive at high temperatures because the Zeeman energy and the energy of the secular part 
of spin-spin interaction are statistically independent in this case. The statistical independence 
denotes that one can change a value of one energy without changing a value of the other, i.e. 
spin system is decoupled to two subsystems the energies of which are additive. 
To explain the latter now we consider a quasi-equilibrium distribution, described by the 
following density matrix: 

Ps = Q7 1 exp{-aH z - pHss] , ?r -, 
Q s = Spexp{-aH z -pU ss } y J> 

Such a form of quasi-equilibrium distribution takes place due to the fact of the availability of 
two invariants of motion. In Equation[25]parameters a and /3 linked to the operators T-i z and 
Hss are thermodynamically conjugative parameters for the Zeeman energy and the energy 
of spin-spin interactions respectively. We can expand the exponent in Equation[25]in powers 
of aH z and and keep only the linear terms. As we shall see later such a linearization 
corresponds to the high temperature approximation. In the linear approximation in aH z and 
pUss, the density matrix is reduced to 

p s = {l-ocn z -pH ss }/Spi (26) 

Here 1 is the unit matrix. Therefore these expressions describe the expectation value of the 
Zeeman energy and the expectation value of the energy of spin-spin interactions: 



2 /c -1 (27) 



(H z ) = Sp P sU z ~ -aSpH 2 z /Spi 
(-Hss) = SppsHss -pSpHl/Spl 

where we have taken into account the orthogonality of operators H z and H ss following from 
Equations |23l and 1241 

SpUtHss = 0. (28) 
One can see from Equations l27l that in the linear approximation in ocH z and pH ss ' 

On the other hand we have the explicit expression for this partial derivative : 

d{H z ) _ d Spexpi-uHg-pHnyHz 



Spexp{-a?T; z — [SHss} 



= -{H z Hss) + {H z ){n ss ) (30) 
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Comparing the right-hand sides of Equation [29] and Eqquation [30] we conclude that in the 
linear approximation in ctH. z and frhL ss 

(U Z U SS ) * (U z ){Hss) (31) 

The Equation [31] means the statistical independence of the first integrals. We remind the 
reader that two or more events are statistically independent if each individual event is not 
influenced by the occurrence of any other and Equation [31] corresponds to the well-known 
rule of probabilities multiplying. 

An entropy of the quasi-equilibrium distribution defined by Equation[26]is equal in the linear 
approximation in olH z and fiH S s to 



S = -k B Spp s lnp s ~ -k B {-]nSpl + u. 2 SpHl/ Spl + p 2 SpHi/Spl}. 



(32) 



We see that the entropy is just a sum of contributions coming from the Zeeman interaction 
and spin-spin interactions. Each of these contributions depends on a or /5 only. The obvious 
identities follow from Equation|32l 



dS 



dS 



1 



(33) 



p k B d(H S s) hT ss 
where the parameters T z and T ss have the physical meaning of temperature of the Zeeman 
subsystem and temperature of the reservoir of spin-spin interactions respectively. The linear 
approximation in txH z and pW S s corresponds to the situation when the energy of spin in 
external magnetic field and one in internal fields are small compared to the heat energy. 
Therefore such a linearization corresponds to the high temperature approximation. Further 
according to the tradition we shall call the parameters a and f> the temperatures of subsystems. 
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Fig. 3. The distribution of the spins over the energy levels of the spin system at high 
temperatures. It is possible to say about two spin temperatures - a and f> in high-temperature 
approximation 

In the thermal equilibrium with lattice the spin system is characterized by one single 
temperature (the temperatures ct and f> coincide). In the general case any thermodynamic state 
of spin system in high temperature approximation can be described by two temperatures: 
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one for a Boltzmann distribution between the Zeeman levels, and a second one describing 
the ordering of the spins in the local fields induced by spin-spin interactions. A pictorial 
illustration is given in Figure|3] 

So at high temperatures we have two statistically independent subsystems - the Zeeman 
subsystem and the reservoir of spin-spin interactions - which are characterized by the 
corresponding independent temp eratures. It i s the approxim ation in which the well-known 
Provotorov's theory was built {Provotorovl Il962allbl: Il963h .. The developed two spin 
temperature formalism approach was appeared to be very fruitful to understand and to 
describe many problems of magnetic resonance, like saturation, sp in-lattice relaxation etc. 
jAbragaml,l2004l;lAtsarkin & RodaUll972l;lGoldmarlll970l;IWolj,ll979h . 

3.3 Low-temperature thermodynamics of spin system 

At low temperatures the factorization condition ( see Equation I3U is violated because in the 
expanding of exponent one has to take high-order terms into account and the right-hand side 
of Equation [29] will not be equal to zero. Consequently ou r subsystems became statistically 
depe ndent ( for more details the reader should refer to jPhilippotl, [1963; iTayurskiil Il989l; 
Il990h ). The advantage of the above-mentioned choice of thermodynamic coordinates is lost. 
Besides the Zeeman temperature doesn't have the physical meaning of temperature, and the 
temper ature of the rese rvoir of spin-spin interactions is the real temperature of the whole 
system jPMlippod.ll%4r) . 

In order to be convinced of it we now note the following. As long as the Zeeman energy 
and the energy of spin-spin interactions are the first integrals then any linear combination of 
these energies with constant coefficients will be the first integral as well. In particular the total 
energy of spin system 

E s = Spp s Hs = Spp s H z + Spp s Hss = (H zz ) + (H ss ) (34) 
and the expectation value of z-component of the total spin 

{s) = spps£s^ = -^spp s n z = -^(n z } (35) 

are the first integrals. If we choose these invariants as thermodynamic coordinates we should 
have an analogy of a thermodynamic system with fixed number of particles. The total number 
of particles and the z-component of total spin are the operators with discrete spectrum and 
these operators play a similar role. Therefore we can use a method analogous to that of the 
grand canonical ensemble and consider a Gibbs ensemble of systems characterized by the 
density matrix 

ps = Q7 1 exp{-/3ft s + vS z } 

Q s =Spexp{-pns+vS z }, w °> 

where the multiplier v is fixed by the condition that (S z ) is given ( now (...) denotes the 

averaging with the density matrix |36l). 

The entropy of the distribution in Eq.[19] may be written as 



S = -k B Spp s \np s = -k B {-)nQ- p{H s ) +v(S z }}. 
Then taking into account the identity 



(37) 



12 



Will-be-set-by-IN-TECH 



(H s }=tu (S z } + {Hss), 
which follows from the Equations[22]and[23]the entropy may be written as 



(38) 



S = k B p(Hss) - k B (S z > (v - wtf) +k B \nQ (39) 

One can use the definition of temperature for systems with fixed number of particles and 
obtain from E quations l37l and 1391 



(S z > 



d{H ss ) 



(40) 



This means that the temperature of the reservoir of spin-spin interactions is the real 
temperature of spin system (see Figure H). On the other hand one can write from Equation 
[37]the expression for the total energy 




Fig. 4. The distribution of the spins over the energy levels of the spin system at low 
temperatures. Only the temperature of the reservoir of spin-spin interactions /3 is the real 
temperature of spin system 



1 n V 



InQ 



whence a simple relation between multiplier v and the chemical potential follows 



(41) 



(42) 



(here we used the definition of the chemical potential for the grand canonical ensemble ). The 
state of the Zeeman subsystem is determined by the value of z-component of total spin or by 
the multiplier v related to the chemical potential. Comparing the Equations |25l l36l and l42l one 
gets a simple relation between the chemical potential of spin system and the parameters a. and 
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For further details conc erning the using of chemical potential in spin thermodynamics the 
reader should refer to jPhirippotl, Il964h . To calculate the expectation values of energies 
("Hz) and (H ss ) at low temperatures one has to take into account the higher-order terms 
in the expansion of the density matrix in Equation [26] As a consequence the factorization 
condition[28]is violated and the Zeeman subsystem and the reservoir of spin-spin interactions 
cannot be considered as independent. So the advantage of the above-mentioned choice of 
thermodynamic coordinates is lost. Besides at low temperatures the entropy written in terms 
a and /5 

S = tx(H z ) + p(Hss) + k B In Q s (44) 

is not a sum of items each of which depends on a or jS only. Therefore at low temperatures the 
parameter a has no physical meaning of the inverse temperature of the Zeeman subsystem. 
In such a situation it is more plausible to use the expectation value of jz-component of total 
spin ( in other words the magnetization of spin system or the number of flipped spins) and 
the total energy as thermodynamic coordinates. Chemical potential and spin temperature 
are thermodynamically conjugated coordinates for them respectively. Such a choice of 
thermodynam i c coordinates ( with taking into account the Equation [43] ) has been made in 
Ide Haas et ailll98d;[l98ll;lTayurskiil [l989; 1990). We note that one can also choose a or /? as 
the thermodynamic coordinates at low temperatures but it is necessary to take into account 
the correlation between them and the fact that parameter a has no physical meaning of the 
inverse temperature of the Zeeman subsystem. 

So, by use very simple, model system - systems of interacting spins in an external magnetic 
field - we have demonstrated that the definition of temperature depends crucially on the 
assumptions made for this procedure. And one has to be very careful doing so. 

4. Conclusion 

The physical quantity "temperature" is a cornerstone of thermodynamics and statistical 
physics. In the present paper the short introduction to the classical concept of temperature 
for macroscopic equilibrium systems was given. The concept of temperature was discussed 
regarding the nanoscale physics and non-extensive systems. It was shown forget that 
it is necessary to remember about the conditions to be satisfied in order to introduce 
"temperature" in macroscopic physics. The concept of "spin temperature" in condensed 
matter physics was reviewed and the advantages of thermodynamic approach to the problems 
of magnetism was illustrated. Partially, two temperatures spin thermodynamics was analyzed 
and the conditions when such approach is valid was studied. 
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